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Abstract 

Integrable multistate or multiflavor/color models were recently introduced. They are gen- 
eralizations of models corresponding to the defining representations of the U q {sl(m)) quan- 
tum algebras. Here I show that a similar generalization is possible for all higher dimensional 
representations. The i?-matrices and the Hamiltonians of these models are constructed by fu- 
sion. The sl{2) case is treated in some detail and the spin-0 and spin-1 matrices are obtained 
in explicit forms. This provides in particular a generalization of the Fateev-Zamolodchikov 
Hamiltonian. 
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1 Introduction 



New one-dimensional integrable lattice models were introduced in Q, within the framework of 
the Quantum Inverse Scattering Method |2], |3], Q . They correspond to a generalization whereby 
every state of the original model, for the defining representation of the Lie algebra A m _\ (see 
for instance ||), is replaced by an arbitrary number of copies. The structure of the #-matrices 
is left unchanged by the replacement, but a usual property, crossing unitarity, is lost. It is 
however still possible to construct integrable open boundary conditions. The eigenvalue set of 
the transfer matrix is unchanged, while the degeneracies increase ||. For integrable periodic 
boundary conditions the eigenvalues and degeneracies change Q . 

These models are of interest in some one-dimensional reaction-diffusion processes @. (Similar 
but different models were also considered in j|] as generalizations of the t-J model.) They also 
appear in connection with the Hubbard model for electrons, as a specific limit (9|, or, for 
m = 2 and at a certain value of the quantum parameter, as building blocks for the natural 
generalizations of the Hubbard model [|ll], 12 1. Upon fermionization, the multiple-states appear 
as different flavors of fermions. The general Hubbard models can then be seen as multichannel 
versions of the original model. 

All the 'multiplicity' A m -\ models of [|l[] correspond to the defining representations of the Lie 
algebras A m -\. As no quantum group formalism H is yet known for the multiplicity models, it 
is not clear whether generalizations to higher dimensional representations are possible. There 
is no obvious or systematic way to do a replacement of states by multiple copies in higher 
dimensional #-matrices. The fusion method, known to work for ordinary models, turns out to 
give the correct answer for the multistate models. 

The main result of this work is to show that such generalizations exist. I first review the 
fusion method for obtaining higher dimensional solutions of the Yang-Baxter equation starting 
from an arbitrary solution. The results of section |2| are quite general and require only a minimal 
number of assumptions. This method is then shown to work for the models at hand. This yields 
multiplicity models corresponding to some higher dimensional representations, and also shows 
that iterated fusions are possible. Thus all higher dimensional representations of sl(m) have 
multistate generalizations. Hamiltonians are obtained as the derivative of the new #-matrices. 
I consider in detail the spin-0 and spin-1 models of sl(2), and give explicit expressions for their 
#-matrices. 



2 Fusion 



The multiple fusion procedure for sl(m) #-matrices was developed in [13]. The idea for other 



algebras is the same 1 14 , 15 ] . A simple description was given in [16] for two successive fusions. 
Here I review the fusion method and give some additional general results. 

The models are defined through their #-matrices. These are solutions of the Yang-Baxter 
equation (YBE): 



#12 (A - (i) #13 (A) # 23 0«) = #23 (M) #13 (A) #12 (A - fi) 



:i 



Here and below, the notation Oij (i ^ j) means that the operator O acts non-trivially on the i 



th 



and j th spaces, and as the identity on the other spaces: Oij = 



(if i < j), where O = J2k a k®bk- Note also that the three spaces 1, 2 and 3 need not be copies 
of the same space. 

Consider any solution # of the Yang-Baxter equation (|J), which becomes proportional to a 
projector at some special value p of the spectral parameter. Define the projector ir^ so that 
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7i"W oc R(p), and let tt^ = I— 7i"W be the orthogonal complementary projector. Setting X—p = p 
in the YBE one obtains: 

7$ R 13 (X) R 23 (A - p) = R 23 (A - p) fl 13 (A) 7$ (2) 

Left and right multiplication of (|2|) by 7r^ yields two equations 

tt^ R 13 (X)R 23 (X-p)7r ( S =0 (3) 
4 2 2 ) i?23(A-p)i?i 3 (A)7ri 1 2 ) =0 (4) 

Define two fused matrices by 

R%>sW = ^R 13 (X)R 23 (X-p)tt$ (5) 
R%>sW = 7T ( §R 13 (X)R 23 (X-p)7T ( g (6) 

Using equations (||,||) and the YBE one shows that the matrices (|5||6|) satisfy a YBE where one 
space is a tensor product of two spaces: 

R% >3 (X-p)R% >4 (X)R 3i (p) = R34(p)R% >i (X)R% >3 (X-p) , i = l,2 (7) 

Thus starting with a given solution of the YBE, we have obtained higher-dimensional i?-matrices 
which are also solutions of the YBE. If di, i = 1, 2, 3 are the dimensions of the spaces 1, 2 and 3, 
then -R<i2>3(A) is a d\ x d 2 X d 3 dimensional square matrix. However the projection operators can 

be diagonalized simultaneously with a change of basis matrix S. 5 , {^ 1 i?^ 12>3 (A) S\ 2 also satisfy 
equation ([?]), and their matrix expressions now contain a number of rows and columns with 
only vanishing elements, and which can be deleted without spoiling the Yang-Baxter property. 
The remaining rows and columns can eventually be relabeled according to the states of the 
corresponding representations. Let tr^ 1 )) = d. Deleting the vanishing rows and columns from 
Si 2 R<\ 2>3 (X) S12 and R<\ 2>3 (X) S\ 2 , yields a d x d 3 and (di x d 2 — d) x d 3 dimensional 
square matrix, respectively. 

It is possible to fuse two matrices -R<i2>3(A) to obtain a matrix -R<i2><34>(A)- Setting p = p 
in (|7|) and multiplying by the projection operators yields two equations similar to (||,||). This 
leads to the following definitions: 

^<12><34>( A ) = ^34 - R <12>4( A + P)- R <12>3( A ) 7r 34 ( 8 ) 
4i2><34>(A) = 4 2 4 ) 4 2 i2>4(A + p)4l2>3(A)4 2 4 ) (9) 

Using the same methods as above one easily shows that these matrices satisfy two Yang-Baxter 
equations (i = 1,2): 

- R <12X34>(A - P) -R<12>5(A) -R<34>5(M) 
= - R <34>5(^) R <12>$W - R <12X34>(A ~ P) (10) 
-^<12><34>(A — p) ^<12><56>(A) ^<34x56>(a0 
= ~X34X56> (M) ^<12><56> W ^<12X34> ~ P) (H) 

Again S^l R^ l2><3A> (X) £34, i = 1, 2, still satisfy the above YBE's. Removing by hand a certain 
number of vanishing columns and rows, one obtains matrices with smaller dimensions. 
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Assume now that the original .R-matrix is regular and unitary, i.e. 



R 12 (0)=cV 12 , i2i2(A)i?2i(-A) = /(A)I , (12) 

where V is the permutation operator. The function /(A) is then even, and c is an arbitrary 
non-vanishing complex number. The fused matrices (||,||) inherit this property. It is however 
necessary to correctly normalize them for the corresponding value of c not to vanish. This 
is achieved by the following normalization: redefine -R<i2>< 3 4>(A) as the right-hand side of 
fll) multiplied by (/(A + p)) . The reason for this normalization is the vanishing of f(p), for 
otherwise the projector R\ 2 {p) would be invertible and therefore equal to the trivial identity 
projector I. 

The limit A —> for the redefined matrices (|8||9|) has to be taken with care. Using the YBE 
one arrives at: 

4 1 L<34>(0) = -c 2 V 13 T2aS^\Ss 12 S^w^Su (13) 

4 2 12><34>(0) = +C 2 V 13 V2AS{ 2 1 7T? 2 ) S 12 S^7T^S U (14) 

^<i2><34>(A)i4 4 ><i2>(-A) = [/(A)] 2 5r 2 1 4 ) 5 12 5 34 1 4 ) 5 3 4 , i = l,2 (15) 

where -R< 34 ><i2>(A) = ^13^24 -R<i 2 >< 3 4> (A) ^13^24- 

Having obtained regular solutions of the YBE, the usual procedure for constructing integrable 
spin-chain Hamiltonians with short-range interaction is by taking logarithmic derivatives of the 
transfer matrices at A = 0. The quadratic hamiltonian density of such integrable hierarchies is 
nothing but the derivative at A = of the matrix R(X) = VR(X). Inclusion of the normalizing 
factor gives: 

£<12><34>(A) = j^ T ^ ) Si2 1 Su 1 4MlRMX + p)Ris(X) (16) 



xi?24(A)i? 23 (A - p)ir{fa®S 12 S M , i = 1, 2 



Taking the limit A — > 0, I find: 



^ f>M f\\ r-r- 2\ / (0) rf-1 (i) a o— l Bo 

^■ jK <12><34>l A )|o - -(TC J 2 J/(o) 12 34 ^ 4 



+ 



1 — i — 1 m m d 2 



2r[Q) SuS£*M ^2 (^32(A + p) J R 13 (A) 



x J R24(A) J R 23 (A-p)) | 7rg7rg5i2^34 (17) 
7 lo 

The two signs are as in the regularity equations above. The first term is proportional to the 
identity and may be dropped. 

It is in fact possible to fuse an arbitrary product of R matrices to obtain solutions of the 
YBE corresponding to most representations of a given Lie algebra. This was carried out in detail 
for sl{2) and si (3) in |L3|, p^| , and also works for the other algebras |l5| . 

The above fusing scheme is now shown to work for the multiplicity A m ^\ models. 



3 New integrable models 

As seen above, apart from satisfying the Yang-Baxter equation, the projector property is the 
only additional ingredient needed to construct fused matrices. In particular R does not have to 
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correspond to a smallest representation to be able to apply the fusion method. The degeneration 
of the generically invertible P-matrix to a projector, for a certain value of the spectral parameter, 
is not automatic. It is however a quite common property, especially for the P-matrices based on 
Lie algebras. For the models studied here, a quantum group structure does not seem to exist, 
and one has to verify that a projector point exists. 

I now apply this formalism to the models of jl]]. This will yield new matrices which define 
new integrable hierarchies. The models of Q are defined as follows. Take positive integers 
and m such that 

m 

2_, n i = n an d 1 < n\ < ... < n m < n — 1 (18) 

8=1 

The inequality restrictions avoid multiple counting of models, but can otherwise be relaxed. The 
set of n basis states is the disjoint union of m sets Ai'. 

card (Ai) = Hi , Ai n Aj = for i / j (19) 

Ai should not be confused with the Lie algebra sl(i + 1). Let E a @ be a square matrix with a 
one at row a and column (3 and zeros otherwise. Define the following operators: 

p(+) = j2 X X Eajai Eaiaj ( 20 ) 

l<i<j<m a.i&Ai otj£Aj 

p(-) = XXX E a ^ ® E a ^ (21) 

l<i<j'<m c«iG-4i (XjdAj 

p( 3 ) = X X X + x~ 1 E a ' a > ®E a w\ (22) 

(23) 

The twist parameters were taken equal to one arbitrary complex parameter x. Let y = e ,A and 
q = e* 7 , where A is the spectral parameter and 7 the quantum parameter. The P-matrix is then 
given by: 

R(\) = V sin(A + 7) + P sin A (24) 
P = P® - (g- 1 P« + g (25) 

This model is denoted by (m, n m ; m, n). For n = m and x = 1 one obtains the An-i 
P-matrix of 0]. P(A) is regular and unitary: 

#12(0) =Pi2 sin 7 (26) 
R12W R21 (-A) = I /(A) = I sin( 7 + A) sin( 7 - A) (27) 

where P2i(A) = Vu Pi2(A) Vi2- As seen in section ||, these properties are inherited by some of 
the fused matrices. 

The right-hand side of (|27| ) shows that A = ±7 are possible projector points. Further checks 
show that matrix ( |24| ) yields a projector at A = p = —7: P(— 7) = — PSU17. Let 

kM = —L^P , vr^I-vrW (28) 

One has: (tt^) 2 = 7rW, i = 1, 2, and vr^ 1 )^ 2 ) = ir^n^ = 0. The dimensions of these projectors 
are given by their traces: 

trfrW) = X^^ = ^(^-XK) 2 ) ( 29 ) 

tr (7r( 2 )) = n 2 -J2mn J = \(n 2 + Y,(m) 2 ) (30) 

i<j \ i / 
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The matrix S which diagonalizes and tt^ is given by: 

m 

S = E Y Y E a ^®E PA 
i=l a,£A fteA 

+ E Y Y (e^ ® E a M - -E a i a i g> £ 



q 



+ E E E (E a > a i®E a ^ + —E a i a *®E a ^\ (31) 

l<i<j<mai£Ai ctj^Aj ^ 9 / 

m 

5-1 = g ^ ^ E a iai ^ E^Pi 

+ — T 2 EE ('-S a * a< <g> - (g) E 

X ~^ X l<i<j<m ai^Ai (XjEAj ^ 

+ — T Y Y Y [E aia i ® E a i ai + —E a i ai ® E aia A (32) 

The diagonalized projectors then read 

S- 1 71-W 5 = E EE ^ w ® - E " ia! ( 33 ) 

s-v 2 )s = E E E 

i=l «,eA fteA 

+ E EE E aiai ®E a ^ (34) 

l<i<j<m cti£Ai otj£Aj 

(i) (i) 

The explicit expressions of the four fused matrices -R< 12 >3 an d ^<i2><34> can be found by 
straightforward if tedious expansions of the products in (|||]J8||]), using the explicit expressions 
(|2^ ) |2l^j3~iy3^) . Similarly, the quadratic hamiltonian density is obtained by replacing the matrices 
in ( |T7|) and expanding the products. 

The following point is worth mentioning. Since Ri2(— 7) -^21(7) = 0, it may seem that the 
complementary projector tt( 2 ) is proportional to ^21(7)- However this is generically not the case. 
It is true only when all the n« are equal to one and x = q . This is another distinguishing 
feature of the nj 7^ 1 models. 



4 Some s/(2) models 



Consider now the m = 2 case, i.e. the XXC models [17]. They have an underlying sl(2) 
structure. Their R matrix is just the multistate version of the one corresponding to the spin-^ 
model (n\ = n<i = 1). 

Tedious but simple calculations lead to the following matrix which carries spin-Ox spin- ^: 
#<12>s(*) = sinA sin(A + 2 7 ) x (35) 
YY X Ea2a2 ® E ^ ai ® EaiPl + E E x^E 13202 ® E aiai ® E a2fi2 

The dimension of this matrix is n\ 712 n, as expected. The corresponding two-dimensional vertex 
model has n\ n<i possible states on, e.g., the horizontal links and n states on the vertical links. 
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Note that the x + x~ l denominator has dropped out of the final result and therefore R is 
defined for all finite values of x, including x = dbi. This also holds for the three fused matrices 
given below. 

The matrix which carries spin-Ox spin-0 can then be obtained: 
4 ><34> (A) = ^(A-7) S in(A + 7)sin(A + 2 7 ) ^ ^E^^E^^E^E^ 

8111(27 ^ «l,A«2,ft 

sin(A - 7) sin(A + 7) sin(A + 2 7 ) ( ^ a) 
~ sin(2 7 - A) /l3 ' 2i [6b) 

and V^ 1 ^ are the permutation operators in the subspaces A2 and Ai, respectively. The 
dimension of this matrix is (ni?i2) 2 , as it should be. It satisfies the regularity and unitarity 
properties (13,15). Here the vertex model has n\ ni possible states on both the horizontal 
and vertical links. The resulting spin-chain is however rather trivial as the A-dependence can be 
normalized away, and the operator part yields S^^u $12 S^tt^ S34, i.e. the identity operators 
in the spin-0 subspaces. 

The matrix carrying spin-lxspin-^ is: 

7T- T^i2>s(A) = +sin(A + 2 7 ) y E a ^ E^ ® E^ (37) 

sm(A + 7) ^ 

+ sin(A + 2 7 ) 5^ E a '^ 2 <g> E~ t2a2 ® E^ 2 

+y sin(2 7 ) ^ ^ E a2ai ® E^ a2 ® E aif32 

"i a 2 ,/3 2 

+qy sin 7 £ £ S^ 1 " 1 ® £° 2/31 ® £ aiC * 2 
+x sin(A + 7) E<3iai ® ^ a2a2 ® 

ai,/3l "2 

sin(A + 7) Yl E ° iai ® jB&Q2 ® Sa2/?2 

"i a 2 ,/32 

+x?T 1 sin7^ 5^ E aia2 ® E 132 ? 2 ® E a2ai 
ai q 2 ,/?2 

+2T W 1 sin(2 7 ) ^ ® J B /3lQ2 (g>£ a2/31 

«2 

+x 2 sin A Y Yl Ea2 ° 2 ® ® 

Q l «2,/32 

+x" 2 sin A Y Yl E ° iai ® ® • BQ2Q2 

ai,/3l "2 

The vertex model model has (ni) 2 + (n2) 2 + ni n2 states on the horizontal links and n states on 
the vertical links. 

The preceding matrix is then used to find the spin-lxspin-1 matrix: 



x ^<12><34>(A) - (38) 



sin(2 7 - A) ^ r(2 ) 
sin(A + 7) 

+ sin(A + 7) sin(A + 2 7 ) ^ £ Ql/31 ® E^ lSl ® <g> £ 5171 



<*i,/3i,7i,5i 
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+ sin(A + 7) sin(A + 2 7 ) ^ E a ^ 2 ® E^ 2 ® E' 32 " 2 ® E 1 ^ 72 

«2,/32, 72,^2 

+y sin(2 7 ) sin(A + 7) ^ i?" 2 " 1 ® E^ 2 ® £ Ql ° 2 ® £T»^ 

«1 «2,/32,72 

+IT 1 sin(2 7 ) sin(A + 7) ^ ^ ® E" 71 " 2 ® E Aai ® E"* 271 

ai,/3i,7l aa 

+y sin(2 7 ) sin(A + 7) ^ ^ E aif31 ® £° 271 ® E' 31 " 1 ® £ 7lC * 2 

ai,/3i,7i Q 2 

+JT 1 sin(2 7 ) sin(A + 7) ^ ^ ^ Ql " 2 ® jB/32 72 ® ® ^ 72/?2 

«1 «2,/32,72 

+ (sin7sin(27) + sin Asin(A + 7)) ^ ^ ^ ai/31 ® ® ® S/32 " 2 

"l:/3l 02,/?2 

+y 2 sin7sin(27) ^ ^ £ Q2Q1 ® ® £ QlQ2 ® E /3l/32 

ai,/?i 02, /?2 

+y" 2 sin7sin(27) ^ ^ E aiQ2 ® E^ 2 ® E a2ai ® E@ 2/31 

a2,/?2 

+2x g _1 y cos 7 sin(27) sin A ^ ^ E" 2 " 1 ® E^ 2 " 2 ® £ ai/31 ® E^ 2 

Ql, /8l «2,/32 

+x- 1 gysin7sinA ^ ^ ® £ a2ai ® E^ 2 ® E^ 2 " 2 

ofl,/3i a2,/?2 

+x 2 sin A sin(A + 7) ^ ^ E aiPl ® E 92 " 2 ® E^ 71 ® £ 7lQl 

ai,/3i,7i "2 

+x 2 sin A sin(A + 7) J2 ® ® E<X ^ ® jE/3272 

Q l "2,/32,72 

sin A sin(A + 7) ^ ^ ® ^ 7iai ® jB/?l71 ® ^ a2a2 

ai,/3i.7i ° 2 

+x- 2 sinAsin(A + 7)^ X! F 1 " 1 ® £ Q2& ® # 72 ® £ 7M 

a l 02,/?2,72 

+x 3 g y" 1 sin 7 sin A ^ ^ E aiC * 2 ® E^ 2 ® £ ft/31 ® £ Q2Ql 

ai>/3i "2,/32 

+2z - V 1 V~ 1 cos 7 sin(2 7 ) sin A ^ ^ E aiai ® E^ 1 " 2 ® E a2/31 ® E^ 2 

«2,/32 

+x 4 sin(A - 7) sin A ^ jE;a2Q2 ® ® 

+x" 4 sin(A - 7) sin A ^ ^ F 1 " 1 ® E^ 1 ® E* 2 " 2 ® £^ 2/32 

This matrix does satisfy the regularity and unitarity properties ( p^Jl5| ). The number of states 
on both the horizontal and vertical links is now {ji\) 2 + (n-2) 2 + n% ri2- 



The four matrices ( 35 , 36 , p7| J3q ) have 2, 1, 10, 19 types of terms, respectively. This was 
expected from the S z conservation of sl(2). When n\ = 712 = 1 the matrices give the si (2) 
2-, 1-, 10-, 19-vertex models. The latter two models can be found in [18, 19, ^Cfl . To carry 
out a comparison it is necessary to relabel the states so that R looks like an operator acting 
on the tensor product of two spaces: R = ^E ® E. For instance, in (f^t), with A\ = {1}, 
A.2 = {2}, one can take the matrix element sin(A + 7) sin(A + 27)£' 11 ® E 11 ® E 11 to correspond 
to I + 1) J + ±) x I + 1) I + 1>, and (x q y)" 1 sin(2 7 ) sin(A + 7)E n ® E 12 ® E 21 to | + 1) | - \) X 1 0> j + \). 
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The matrix (19) of reference |2(| can be obtained by fusing the symmetric form of the sl(2) R- 
matrix fj~7|| . It is also necessary to do a gauge transformation after fusing, in order to render 
the resulting matrix completely symmetric. This introduces the unusual square-root element d 
of that reference. Finally, let u = A + ry/2 and r/ = 7 to complete the identification. Similar 
remarks apply to the identification of the 19-vertex models. 

The interpretation of the multiple-states in terms of s/(2) states is simple. Recall first that 
for the models (|24|) with m = 2, the interpretation is: 

state(ai) < — ► | + -) ai , state(a 2 ) < — ► I - ^)a 2 ( 39 ) 

This yields the following identifications when two spin-^ spaces are fused: 

state(a 2 ai) < — > |0) a2ai , nira 2 (40) 

for the spin representation obtained from n^, and 

state(«i/3i) < — > I + l)ai/9i , {nij 2 (41) 
state («i a 2 ) < — ► |0} ai a 2 , ni ra 2 (42) 
state(a 2 /3 2 ) < — ► | — l) a2 /3 2 , (n 2 ) 2 (43) 

for the spin 1 representation obtained from tt^ 2 \ The numbers on the right are the number 
of states of the corresponding type. Whereas the numbers of copies of the states | ± 1) are 
uncorrelated, (ni) 2 and (n 2 ) 2 , the number of states of type |0) has to be ri\n%. Thus a naive 
trial to obtain the multistate version of the spin-lxspin-1 model would have failed. Fusion gives 
the correct answer. 

The first derivative at A = of -R<i 2>< 34>(A) = ^13^24 -R<i 2 ><34>(^)> for the matrix of 
(|38|), gives the hamiltonian density Ha + i for the multistate version of the spin-1 model. The 
n\ = i%2 = 1 Hamiltonian is also known as the Fateev-Zamolodchikov model Ipi 



5 Concluding remarks 

The multiplicity A m _i models were shown to allow fusion and multistate models were obtained 
for higher dimensional representations. The spin-0 and spin-1 models were derived explicitly. 
This provided in particular the generalization of the Fateev-Zamolodchikov model. 

The underlying sl(m) structure of these models was justified in It provides a natural 
way to label them in terms of sl{m) representations. This language has been used throughout 
the paper and in particular in section ||. I recall here the main points in view of an extension 
to the higher dimensional models just obtained. The matrices ( ^4[) have the following structure: 
R{\) = J2k9k{ty®ki where all the spectral parameter dependence is contained in the functions 
gk(X). The operators Ok may depend on the other parameters. Replacing R in the YBE's 
it is required to satisfy, and identifying the coefficients of the linearly independent functions 
yield a set of equations for the operators Ok, where the spectral parameters do not appear. In 
the case of the Ok of (|24| ) one has simple raj-independent trigonometric functions, and obtains 
the Hecke algebra, for all values of the ra^'s, and in particular for m = 1. This and the fact 
that the operators have the same structure as for = 1 shows that one has the natural 
generalization. The algebraic Bethe Ansatz was also found to be based on the Dynkin diagram 
of A m -i. Thus the functional and operatorial structures of R(X) = J2k 9k(X)Ok are the same as 
for rii = 1. The quite general fusion construction reviewed in section |2| preserves the functional 
and algebraic structures. Therefore one again has the natural multistate generalization for the 
higher dimensional representations. 
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The rational limit of all the matrices considered above exist and provide rational solutions 
to the various Yang-Baxter equations. This limit is obtained by rescaling A to 7A, and dividing 
by an appropriate power of 7 or sin 7 in order to obtain a finite limit. 

The models (§4|) have extended symmetries |J: sl(ni) © ... © sl(n m ) © u(l) © ... © u(l). It 
is obvious that the higher dimensional models inherit similar symmetries. These are related to 
the existence of many states of a same type, as seen in the spin-0 and spin-1 examples. Most 
(but not all) diagonal operators commute with the i?-matrix, as they either correspond to k(1) 
charges of states of the same type, or to realizations of S z . The non-diagonal operators which 
commute with the i?-matrix are those which acts within the space spanned by a state and its 
copies. 

These symmetries will also be reflected in a diagonalization by Bethe Ansatz, as happened 
for the fundamental models [fflj. The diagonalization makes use of the standard techniques 
associated with higher dimensional representations. In particular, it uses fusion to find relations 
between the various transfer matrices. 
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